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On Yukawa quasi-unification with /x < 
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Although recent data on the muon anomalous magnetic moment strongly disfavor the constrained 
minimal supersymmetric standard model with /i < 0, they cannot exclude it because of theoretical 
ambiguities. We consider this model supplemented by a Yukawa quasi-unification condition which 
allows an acceptable 6-quark mass. We find that the cosmological upper bound on the lightest 
sparticle relic abundance is incompatible with the data on the branching ratio of b — > sj, which is 
evaluated by including all the next-to-leading order corrections. Thus, this scheme is not viable. 

PACS numbers: 12.10.Kt, 12.60.Jv, 95.35.+d 



The constrained minimal supersymmetric standard 
model (CMSSM), which is a predictive version of the 
minimal supersymmetric standard model (MSSM) based 
on universal boundary conditions, can be further re- 
stricted by being embedded in a supersymmetric (SUSY) 
grand unified theory (GUT) with a gauge group contain- 
ing SU(4)c and 577(2) R . This can lead [jj to Yukawa uni- 
fication |2j, i.e., the exact unification of the third gener- 
ation 'asymptotic' Yukawa couplings. However, for both 
signs of the MSSM parameter fi, the CMSSM supple- 
mented by the assumption of exact Yukawa unification 
yields unacceptable values of the 6-quark mass. This is 
due to the generation of sizeable SUSY corrections to 
mb which have the sign of /i. The predicted tree- 
level rrib(Mz), which turns out to be close to the upper 
edge of its 95% confidence level (c.l.) experimental range 
2.684 - 3.092 GeV, receives, for fi > (fi < 0), large 
positive (negative) corrections which drive it well above 
(a little below) the allowed range. This range is derived 
from the 95% c.l. range for mb(rrib) in the MS scheme 
3.95 — 4.55 GeV J3 evolved up to Mz in accord with the 
analysis of Ref. yj and with a s (Mz) = 0.1185. 

In Ref. 0, we have built SUSY GUT models based on 
the Pati-Salam gauge group SU(4) C x SU(2) L x SU(2) R 
which moderately violate Yukawa unification and, thus, 
can allow an acceptable 6-quark mass for both signs of 
fi even with universal boundary conditions. These mod- 
els yield a class of 'asymptotic' Yukawa quasi-unification 
conditions (YQUCs) replacing exact Yukawa unification. 
Supplementing the CMSSM with one of these YQUCs, 
we achieved 6], for fi > 0, an acceptable rrib(Mz) in ac- 
cord with all the phenomenological and cosmological re- 
quirements. The corresponding SUSY GUT model leads 
to a new version of shifted hybrid inflation 

Here, we will examine whether the CMSSM with a 
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YQUC can yield an acceptable rrib(Mz) in accord with 
all the available constraints for \i < too. In this case, 
compatibility of the 6-quark mass with the data requires a 
deviation from Yukawa unification which is much smaller 
than the one needed for fi > 0. Consequently, there is no 
need to extend the simplest SUSY GUT model of Ref. 0, 
which yields a suppressed violation of Yukawa unification 
originating from non-renormalizable interactions. 

The CMSSM with /x < is now strongly disfavored 
(but not excluded (9j) by the constraint on the deviation, 
5a M , of the measured value of the muon anomalous mag- 
netic moment a M from its predicted value in the standard 
model (SM) g™. The recent data [13 on a u confirm the 
earlier ones [ll| but their precision is twice that of the 
previous data. On the other hand, a^ M is not yet sta- 
bilized. This is mainly due to the hadronic vacuum po- 
larization contribution. According to the most updated 
|12| evaluation of this contribution, there is a consider- 
able discrepancy between the results based on e+e~ data 
and on r data. Taking into account both results, we ob- 
tain the 95% c.l. range -4.7 x 10~ 10 < Sa^ < 56 x 10~ 10 , 
where the negative values are allowed only by the r-based 
calculation. Using the formulae of Ref. [13| , we calculate 
8a^ in the CMSSM and we observe that the result is nega- 
tive (positive) for \x < (// > 0) and increases (decreases) 
as the lightest sparticle (LSP) mass tolsp increases. Con- 
sequently, a lower bound on tolsp can be derived from 
the lower (upper) experimental limit on Sa^ for \i < 
(/i > 0). Also, we see that \i < is allowed only by the 
r-based calculation. 

The CMSSM with fi < is, also, severely restricted 
by the recent experimental results |14| on the inclusive 
branching ratio BR(6 — > s-f). We compute this branching 
ratio by using an updated version of the relevant calcu- 
lation in the micrOMEGAs code Comparing care- 
fully this version with our calculation of BR (6 — » .57) in 
Ref. , we concluded that this improved version is com- 
plete and reliable. The SM contribution is calculated 
by using the formalism of Ref. [16j and incorporating 
the improvements of Ref. 01 • The charged Higgs boson 
{H^) contribution is evaluated by including the next- to- 
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leading order (NLO) QCD corrections from Ref. [Tm and 
the tan/3 enhanced contributions from Ref. [lSj. The 
dominant SUSY contribution includes resummed NLO 
SUSY QCD corrections from Ref. which hold for 

large tan/3. The contribution interferes construc- 
tively with the SM contribution, while the SUSY one in- 
terferes constructively (destructively) with the other two 
contributions when fx < (/i > 0). The SM plus H ± and 
the SUSY contributions decrease as ttzlsp increases and 
so, for n < 0, an additional lower bound on tolsp can 
be derived |2(J] from the upper experimental limit on this 
branching ratio. At 95% c.L, the experimental range of 
BR(6 — > S7) is estimated 6] by combining appropriately 
the various experimental and theoretical errors involved. 
The result is 1.9 x 1(T 4 < BR(6 -> sj) < 4.6 x 1(T 4 . 
We find that the NLO corrections reduce BR(6 — > sj) by 
10 — 15%. In Ref. 21] (which adopts the opposite sign 
convention for fi), where the tan/3 enhanced and NLO 
SUSY QCD corrections were ignored, the reduction was 
considerably higher yielding a less stringent restriction. 

The LSP mass is bounded above by the requirement 
that the LSP relic abundance SIlsp^ 2 in the universe 
does not exceed the upper limit on the cold dark matter 
(CDM) abundance derived from observations. Recent re- 
sults [22 suggest that, at 95% c.L, n LSP h 2 < 0.22. In 
general, Olsp/i 2 increases with tolsp- So, its value is ex- 
pected to be quite large due to the large tolsp's required 
by the constraints from Sa^ and, especially, b — > S7. The 
^Lspfr 2 is calculated by using micrOMEGAs [l^, which is 
the most complete code available at present. It includes 
accurately thermally averaged exact tree-level cross sec- 
tions of all possible (co) annihilation processes and treats 
polo effects properly with one-loop QCD corrected decay 
widths. In the model which we will consider, the LSP 
is an almost pure bino (x) and the next-to-LSP (NLSP) 
is the lightest stau (t 2 ). Important annihilation chan- 
nels are not only the ones with fermions // in the fi- 
nal state, but also the ones with HZ, W ± H T , hA j2^| 
(h, H are the CP-even neutral Higgs bosons and A the 
CP-odd neutral Higgs boson). On the other hand, as we 
can observe from the sparticle spectrum of the model 
(see Fig. QJ, is much smaller than 2tolsp- So, the 
LSP annihilation via the s-channel exchange of an A- 
boson is not enhanced as in the \i > case [(J |24| . As a 
consequence, the only available mechanism for reducing 
^lsp^ 2 is bino-stau coannihilation |rl I2HI25I l2r| (for an 
updated analysis, see Ref. [27]]), which is activated when 
the relative mass splitting between the T2 and the LSP, 
Af 2 = (mf 2 — tolsp )/'7iLSP) is reduced below about 0.25. 
For fixed Tolsp, ^lsp^ 2 decreases with A f2 , since coan- 
nihilation becomes more efficient. So the CDM bound 
yields an upper limit on Af 2 . 

We will consider here the CMSSM with fi < which 
results from the simplest SUSY GUT model constructed 
in Ref. |(J without the extra Higgs superfields </>, </>, <//, 
(j) 1 . To reduce the number of free parameters, we will fur- 
ther assume that the non-renormalizable SU(2)r triplet 
coupling of Ref. 6] which leads to violation of Yukawa 
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FIG. 1: The mass parameters raj, mo, M 1 / 2 and Msusy 
versus m L sp for A = 0, Af 2 = and m b (M z ) = 2.888 GeV. 

unification is suppressed. The violation is then achieved 
predominantly via the non-renormalizable SU(2)r sin- 
glet coupling of Ref. [a] and, thus, a\ = 012 in Eq. (15) of 
this reference, with a\, of order 0.1 or smaller. Note 
that the other alternative with the triplet coupling dom- 
inating, which was studied in Ref. with fi > Q, is not 
viable for ji < since it leads to a charged LSP. Under 
the above assumptions, the YQUC takes the form: 

h t :h b :h T = \l-c\:\l-c\:\l + 3c\, (1) 

where c is a complex parameter of order 0.1 or smaller. 
We will take c < 0, where the maximal splitting between 
/if,, h T and, thus, the largest enhancement of rrib(Mz) 
is achieved for fixed |c| 0. For —1/3 < c < 0, this 
splitting is — 5h = (hb — h T )/hb — — 4c/(l — c). We see 
that the tau Yukawa coupling becomes smaller than the 
top and bottom Yukawa couplings which remain unified, 
in contrast to the /i > case of Ref. 0, where the tau 
and bottom Yukawa couplings split from the top Yukawa 
coupling by the same amount. The SUSY GUT model 
used leads |8] to successful shifted hybrid inflation. 

Below the GUT scale M GU t, this SUSY GUT model 
reduces to the CMSSM with universal soft SUSY break- 
ing scalar masses too, gaugino masses M1/2, and trilin- 
ear scalar couplings Ao at Mqut- It also implies the 
YQUC in Eq. We closely follow the renormaliza- 

tion group and radiative electroweak breaking analysis 
of Ref. ;6| . The conditions for electroweak breaking are 
imposed on an optimized scale Msusy = ( m t 1 m t 2 ) 1 ^ 2 
(ti t 2 are the stop mass eigenstates) , where the spar- 
ticle spectrum is also evaluated. We incorporate the 
two-loop corrections to the CP-even neutral Higgs bo- 
son mass matrix by using FeynHiggsFast [28j and the 
SUSY corrections to TOt, (Msusy) and m T (Msusy) from 
Ref. [2^. The corrections to to t (Msusy) lead |2l| to 
a small enhancement of tan/3. For m t (m t ) — 166 GeV, 
m T (M z ) = 1.746 GeV, and any given m b (M z ) in its 95% 
c.L range (2.684 — 3.092 GeV), we can determine the pa- 
rameters c and tan /3 at Msusy so that the YQUC in 
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FIG. 2: Restrictions on the ttilsp — Af 2 plane for Aq = and 
2.684 GeV < m b (M z ) < 3.092 GeV. The solid (dot-dashed) 
line is the lower bound on ttilsp from 5a M > —4.7 x 10 -10 
(BR(6 -» s-y) < 4.6 x 10~ 4 ). The dashed (double dot-dashed) 
line is the upper bound on ttilsp (Af 2 ) from tt-LSph 2 < 0.22. 



FIG. 3: Restrictions on the ttilsp — Ao/M\/2 plane for any 
Afj, and 2.684 GeV < m b {M z ) < 3.092 GeV. The notation is 
generally as in Fig. [5] The dotted line and the inclined parts 
of the solid line are explained in the main text. 



Eq. (JIJ is satisfied. We are then left with only three free 
input parameters mo, M 1 / 2 and Aq. The first two can 
be replaced [2lJ by ttilsp and Af 2 . In Fig. ^ we present 
the values of rriA, ttto, M l / 2 and Msusy versus ttilsp for 
A = 0, A f2 = 0, m b {M z ) = 2.888 GeV. These values 
are affected very little by varying m b (Mz). Our scheme, 
unlike the /i > model of Ref. Q , yields 2tttlsp 3* m A ■ 

The restrictions on the tttlsp — Af 2 plane for Aq = 
and any m b (Mz) in its 95% c.l. range are given in Fig. EI 
The lower bound on tttlsp from 5a^ > —4.7 x 10~ 10 
corresponds to m b (Mz) — 3.092 GeV and is repre- 
sented by a solid line. The maximal Af 2 on this line 
is about 0.0105 and yields ^lsp^ 2 — 0.22 for the same 
value of m b (Mz). As m b (Mz) decreases, the maximal 
Af 2 from Olsp^ 2 ^ 0-22 increases along the double 
dot-dashed line and reaches its overall maximal value 
A f2 ~ 0.0142 at m b (M z ) ^ 2.684 GeV. The upper 
bound on tttlsp from SIlsp^ 2 0.22 is achieved at 
m b (Mz) — 2.684 GeV and corresponds to the dashed 
line. In the lightly shaded region allowed by 5a^ and 
CDM considerations, 612.4 GeV < tttlsp ^ 873.4 GeV. 
The maximal tttlsp is achieved at Af 2 = yielding 
BR(6 — > s"f) ~ 5.26 x 10~ 4 . The lower bound on tttlsp 
(dot-dashed line) from BR(6 — > sj) < 4.6 x 10~ 4 cor- 
responds to m b (Mz) — 3.092 GeV. In the correspond- 
ing allowed (dark shaded) area, tttlsp ^ 1305.04 GeV 
with the minimal tttlsp achieved at Af 2 = and yielding 
^Lsph 2 — 0.65. So, for A = 0, there is no region where 
all the restrictions are satisfied. Note that the constraints 
BR(6 -> s 7 ) > 1.9 x 10~ 4 and 5a M < 56 x 10~ 10 always 
hold for the CMSSM with /j < 0. Also, m h > 114.4 GeV 
[30| is valid due to the heavy SUSY spectrum. 

Departure from the Aq = case does not change our 
conclusion as we can see from Fig. [3J where the restric- 
tions on the tttlsp — Aq/M-l/2 plane are presented for 
any Af 2 and any allowed m b (Mz)- Actually, we can just 



take Af 2 = 0, where all the constraints happen to be- 
come as less restrictive as possible for any given Aq (see 
Fig. [2J • We follow the same notation for the various lines 
and areas as in Fig. 2. For -1.84 < A /M 1/2 < 1.21 
(-1.67 < A /M 1/2 < 1.03), the solid (dashed) line 
corresponds to Sa^ ~ -4.7 x 10~ 10 (ft LS p/i 2 ~ 0.22) 
with m b (M z ) ~ 3.092 GeV (2.684 GeV). As \A \ in- 
creases, ttia decreases and reaches its lowest allowed value 
(~ 100 GeV for the tan/3's encountered here j33) on 
the boundaries of the lightly shaded regions for fixed Aq . 
Allowing m b (Mz) to increase from its lower value, we 
then obtain the upper and lower inclined parts of the 
dashed line corresponding to ttt^ ~ 100 GeV. Their end 
points lie at A0/M1/2 ~ 1.21,-1.84, where m b (M z ) is 
maximized. At the maximal (minimal) tttlsp on the dot- 
ted line (inclined parts of the solid line), JIlsp^ 2 < 0.22 
(Scift > —4.7 x 10~ 10 ) is not saturated, due to the com- 
bined requirements ttt.4 > 100 GeV and m-f 2 ^ ttt^. 

The maximal tttlsp in the lightly shaded area of Fig. 
which is allowed by Sa^ and CDM considerations, and for 
A /M 1/2 > (< 0) is tttlsp - 868 GeV (967.3 GeV) and 
is achieved at A /Mi/ 2 ~ 1.03 (—1.67). This corresponds 
to the upper (lower) corner of the dashed line, where 
BR(b -» s-y) ~ 5.58 x 10~ 4 (5.45 x 10~ 4 ). Note that, for 
Ao/Mi/ 2 < 0, the processes f 2 f 2 — > W ± H^ 1 1 ZA become 
more efficient. Their total contribution to the effective 
cross section is 3 — 6.5% as A a /Mi/ 2 decreases from 
to —1.67. So, coannihilation is strengthened and bigger 
tttlsp 's are allowed than for Aq/Mi/ 2 > 0. The minimal 
tttlsp in the dark shaded area, allowed by b — > S7, is 
tttlsp — 1295 GeV and is achieved at A /M 1 / 2 ~ —0.2, 
where r^SP^ 2 — 0.62. Thus, even for A ^ 0, there is no 
range of parameters allowed by all the constraints. 

In the lightly shaded area of Fig. |3 tan (3 (c) ranges 
between about 44.1 (-0.09) and 51.5 (-0.007). These 
values are achieved at the lowest corner of this area 
(A /M 1/2 ~ -1.89, m b (M z ) 3.092 GeV) and the 
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upper corner of the dashed line [Aq/Mx/2 — 1.03, 
in-biMz) — 2.684 GeV) respectively. Consequently, the 
splitting — Sh ranges from about 0.33 to 0.028. However, 
hxing mi,(Mz) to its central value, the range of — Sh is 
reduced to 0.13 — 0.09 corresponding to 47.6 < tan/? < 49 
and -0.038 < c < -0.024. 

The present investigation is an improved version of our 
analysis in Ref. 2JJ (the sign of \i there is opposite to the 
one adopted here). The main improvements are the re- 
placement of Yukawa unification by the YQUC of Eq. 
in connection with the more stringent present bounds 
on mj,, the consideration of the <5a M constraint |l3|. the 
inclusion of the NLO SUSY QCD and tan/3 enhanced 
corrections in BR(6 — > sj), and the evaluation of 
^Lspft 2 by the updated code of Ref. ^J. These improve- 
ments turn out to be of crucial importance overruling 
our previous conclusion. Our results are essentially unaf- 
fected by allowing a s (Mz)to vary in its 95% c.l. range, 
which slightly enlarges ii the range of m^Mz). This 
is due to the fact that, here, the LSP annihilation to bb 



via an ^4-pole exchange is subdominant in fl^sph 2 . 

In summary, we studied the CMSSM with /i < and a 
YQUC from the simplest SUSY GUT of Ref. @. We im- 
posed the constraints from mj, <5a M , b — > sj and CDM. 
Although the NLO corrections to b — > sj (coannihi- 
lations) drastically reduce (enhance) the lower (upper) 
bound on tolsp, the b — > sj and CDM requirements re- 
main incompatible. Thus, despite the fact that, with the 
r-based calculation of a^ M , the 5a ^ and CDM criteria 
can be simultaneously valid, this model is excluded. 
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